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Abstract 

In the framework of the effective field theory for interacting BFKL pomerons, applied to nucleus- 
nucleus scattering, boundary conditions for the classical field equations are discussed. Correspondence 
with the QCD diagrams at the boundary rapidities requires pomeron interaction with the participating 
nuclei to be exponential and non-local. Commonly used 'eikonal' boundary conditions, local and linear 
in fields, follow in the limit of small QCD pomeron-nucleon coupling. Numerical solution of the classical 
field equations, which sum all tree diagrams for central gold-gold scattering, demonstrates that corrected 
boundary conditions lead to substantially different results, as compared to the eikonal conditions studied 
in earlier publications. A breakdown of projectile-target symmetry for particular solutions discovered 
earlier in |13j is found to occur at roughly twice lower rapidity. Most important, due to a high non- 
linearity of the problem, the found asymmetric solutions are not unique but form a family growing in 
number with rapidity. The minimal value for the action turns out to be much lower than with the eikonal 
boundary conditions and saturates at rapidities around 10. 
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1 Introduction 



In the framework of the high-colored QCD, in the Regge kinematics and perturbative domain Aqcd « 
\t\ « s, hadronic collisions are described by the exchange of BFKL pomerons (see the review in [T]) 
interacting via the triple pomeron vertex [21 El SI E] - In applying this mechanism to concrete reactions one 
has to distinguish between purely hadronic and nuclear participants. For nuclei the emitted or absorbed 
pomerons are actually interacting with a spatially widely distributed nucleons, so that the color dipole 
density is small on the average. The interaction with individual nucleons is however enhanced by factor 
A 1 / 3 . As a result we have a strong interaction with a dilute colour dipole source ('the dilute regime' in 
the terminology of [6]). This greatly simplifies the treatment since in allows to neglect the contribution 
from pomeron loops, the fact which was stressed already in [7] in the framework of the old local pomeron 
Regge-Gribov field theory (RGFT). In contrast, in a hadron the color dipoles are concentrated in a small 
volume. In this dense regime the contribution of pomeron loops is not damped at all and in all probability 
has a decisive influence. Up to now there exists no reliable method to study the coupling of the pomerons 
to a dense source of colour dipoles nor the loop contribution (see however attempts in |8j). Therefore main 
positive results have been obtained in application of pomeron dynamics to reaction with nuclei. The most 
consistent one is the study of DIS on nuclei described by a well-known Balitsky-Kovchegov (BK) equation 
[Hi [TO]. With certain reservations this equation can also be applied to hadron-nucleus scattering. 

A natural extension of these applications is to nucleus-nucleus scattering. As compared to the hadron- 
nucleus case, the corresponding amplitude already in the tree approximation contains many more diagrams 
not taken into account in the Balitsky-Kovchegov equation, which only sums pomeronic fan diagrams. The 
equation summing all tree diagrams for nucleus-nucleus scattering was formulated and studied numerically 
in papers [11} I12j . In paper [13] the ideas of [111 [TO] were applied for the scattering of two protons. In 
the diagrammatic way solutions for the pomeronic fields, found in these papers, may be represented as all 
possible tree diagrams of the type shown in Fig. [T^, which are enhanced as compared to loop diagrams 
Fig. db by factors A 1 ^ and B 1 ^ for the colliding nuclei with atomic numbers A and B. The interaction of 
pomerons with the target and projectile nuclei was taken in these papers just as a collection of interactions 
with individual nucleons shown in Fig. [Ik, This 'eikonal'-type interaction borrowed for the Glauber picture 
for hadron-nucleus scattering has been commonly used in the study of reactions with nuclei (see e.g. old 
papers dedicated to the zero-dimension variant of RGFT local pomeron models [14 } 115 } 16j) 

However, as first clearly stated by G.A.Winbow [T7J, the eikonal form of the interaction suitable for 
the description of hadron-nucleus scattering is inadequate for nucleus-nucleus scattering. The reason is that 
the effective number of collisions in the nuclear-nuclear interaction is oc A 4 / 3 and so much higher than in 
the hadron nuclear interaction oc A 1 / 3 . Taking this into account leads to a change in the pomeron- nucleus 
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Figure 1: Examples of diagrams of the effective field theory of QCD pomerons interacting with triple pomeron 
vertices: a) a tree diagram defining the classical limit; b) a diagram with quantum loops. 

interaction. It now becomes non-local on the nuclear scale and not represented by pomeron- nucleon vertexes 
like shown in Fig. [2^ but rather by pomeron-nucleus vertexes as a whole shown in Fig. [2b. 

In this paper we study the nucleus-nucleus interaction with this corrected pomeron-nucleus coupling. 
We also compare the following results with those obtained with the old-fashioned eikonal interaction. We 
do not aim here at calculating the full nucleus-nucleus amplitude at all impact parameters needed for the 
AB cross-section. Our primary task is to see how important is the use of correct boundary conditions 
corresponding to pomeron-nucleus interaction. Therefore we limit ourselves to the study of central collisions 
of identical nuclei. 

The paper is organized as follows. In the next section we describe the formalism of the effective field 
theory of interacting BFKL pomerons. In Sec. 3 initial conditions are discussed which follow from different 
choices of pomeron interaction with the nucleus. In Sec. 4 we present semi-classical solutions of the theory 
with different initial conditions and the corresponding values of the S-matrix. We discuss our results and 
conclude the paper in Sec. 5. 



2 Effective field theory for interacting pomerons 

In the perturbative QCD with a large number of colours the pomerons can be described by two real fields 
4>(y,k,f3) and (jy(y,k,l3) depending on rapidity y, relative transverse momentum k of the two reggeized 
gluons which form the pomeron and transverse point j3. In the nucleus-nucleus scattering problem the 
transverse momentum carried by pomerons in tree diagrams is negligible as compared to the relative gluon 
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Figure 2: Vertices for the pomeron-nucleus interaction: a) single interaction with a nucleon, local on the 
nuclear scale; b)multiple interaction with a nucleon, non-local on the nuclear scale 

momentum inside the pomeron. The action for this simplified case was introduced in [T2"] (see [T8J for 
a general case with loops). It can be written as a sum of three terms: 

A = A +Ai + Ae- (1) 

The free part is 

A = 2£dy J d 2 p(^\K^ + H^j \<f>), (2) 

where H is the BFKL forward Hamiltonian for the so-called semi-amputated amplitudes pQ and K is a 
differential operator in k commuting with H 

K = Vlk^Vl (3) 

Symbol (...) means integrating over k with weight l/(2-7r) 2 . Action Aq generates propagators which are the 
BFKL Green functions with operators K~ l attached at their ends times 5 2 {b\ — 62)- The interaction part 
Ai describes splitting and merging of pomerons: 

Ai = ^[c£ dy J fPp^Kj + cpK^). (4) 

Finally Ae is the external action which describes interaction of the pomerons with the colliding nuclei. 
This part is the main subject of discussion in this paper, to which the next section is devoted. 

The form of the action in terms of the fields <fi and $ is directly related to the diagrammatic picture 
of the pomeron interaction. However it is not convenient for numerical studies due to large number of 
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derivatives in the operator K. An alternative description, more easily tractable numerically, can be made 
in terms of the corresponding unintegrated gluon densities, f(y,k 2 ,/3) and p(y,k 2 ,/3) introduced in [19J: 



f(y,k 2 ,(3) 



N, 



c iAx-r2 



4a s 7T 2 



with the inverse relation 
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The form of the action at a given impact parameter in terms of the unintegrated gluon density was 
presented in [13J. Since it serves as our main calculational tool, we rewrite it here in an explicit form for 
clarity. We present the free and interaction parts of the action in the form of integrals over the rapidity and 
transverse coordinates of the Lagrange functions: 
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In terms of / and f* functions Cq j are given by the following expressions, in which we suppress the 
dependence on /3 entering parameter 
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where /Co is the amputated forward BFKL kernel given by 
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In the interaction term we also suppress the y dependence, which enters as a parameter. 
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From this action we standardly obtain the following classical equations of motion for fields / and f' at 
given y and j3 : 
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and a similar equation for f* which is obtained from Eq. (Ilip by the change d y — > —d y and interchange 
f The functional derivatives over the source terms 5AE/5f^(y, k 2 , (5) and 5AE/5f{y, k 2 ,j3) determine 

initial values of the fields / and p at rapidities y = and y = Y at which the target and projectile nuclei 
are moving. 

The external action Ae depends on the overall impact parameter B of the collision. As a result, the 
total action A is also B dependent. At a given B the scattering matrix for the collision is determined by 
the functional integral: 

f DfDf exv(-A{f,f;Y,B}) 

S(Y,B) = -> <-. 12 

JD/Dp e W (-A {f,P;Y,B}) 

Solving the equations of motions for the fields / and f\ ( (jlip and a similar one for p ), we obtain the 
S-matrix in the semi-classical approximation 

S(Y, B) = exp ( - A{f, / f ; Y, B}) , (13) 

where / and p fields in (|13p are to be understood as the solutions of the equations of motion ( (jlip and the 
one for p). The elastic amplitude for the scattering of the two nuclei in this approximation is given by: 

A el (Y, B) = 1 - exp ( - A{f, /t; Y, B}) . (14) 
3 Boundary conditions for the nuclear scattering problem 

The form of the external action is determined by the actual mechanism of the pomeron interaction with 
nucleons in the nucleus. The simplest assumption one can make is to assume that this interaction is similar 
to the one between the pomeron and virtual photon in DIS. Then the vertex for it will be just a convolution 
of the pomeron wave function with the colour dipole density of the nucleon times the probability to find 
the nucleon at a given transverse point (3 described by the nuclear profile function. For the interaction with 
nucleus A moving at rapidity y = the vertex will be 

-T A (/3) <p^(y = 0,P) > . (15) 

Here p{k) is the dipole density of the nucleon and < ... > as before means integrating over k with weight 
l/(27r) 2 . The minus sign comes in the course of transition from the physical complex field to our real field, 
which involves factor i. Likewise the vertex for the pomeron interaction with a single nucleon from nucleus 
B at rapidity Y and transverse point (3 will be 

-T B (B-P)<p<t>(y = Y,P)> . (16) 
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Now we take into account that from each nucleus different numbers of nucleons located at different transverse 
points may interact with pomerons. As a result the S*-matrix for the collision of two nuclei A and B will be 
written double sum 

S{Y,B) = 

A B . ua n B 

E E <°l / II II rf 2 A^T A (ft)r B (B - ^)(- < My = ^,^) > )(- < = o,ft) > )|o), (17) 

n^=0 ng=0 i=lj=l 

where the pomeron Green function is to be taken with the action which describes only their movement and 
mutual interaction: Aq + Ai. Note that the actual interaction involves only terms with ua,tlb ^ 0. The 
added terms with ua = 0, «b / or = 0, ua 7^ are zero and the term with ha = ub = gives unity, 
which converts the amplitude into the S-matrix. Doing the summation over ua and tib we get 

S(Y,B) = (0\(J d 2 (3T A {b){l- < fxj>\y = 0,/3)) A ( J d 2 (3T B {B - (3){l~ < p(p{y = Y,(3)) B \0). (18) 

This corresponds to the external action 

A E = A In (l - J d 2 (3T A ((3) < ptf(y = 0,(3) > ) + Bhi (l - J d 2 (3T B (B - (3) < p<j>(y = Y,(3)>). (19) 

Note that, in contrast to the standard Glauber formulas, under the sign of logarithm the interaction term 
appears integrated over (3. This makes the action nonlocal on the nuclear scale: expansion of the logarithm 
generates multiple integrations in (3. Also integration over (3 compensates the smallness of Ta b at large A 
and B. So, unlike the hadron-nucleus case in the Glauber picture, further simplification of this expression 
is not possible unless the interaction itself is small. Formally p is of order g 2 and, assuming g to be small, 
in the leading order one has to take only the first non-trivial term of the expansion of the logarithm. This 
brings us to the commonly used 'eikonal' form of the external action 

A e = -aJ d 2 f3T A ((3) < P ^(y = 0,(3) >-B J d 2 (3T B {B - (3) < p^y = Y,(3)> (20) 

and returns us to the earlier studies in [T2"| [15] . 

However the coupling constant which appears in p is in fact not the same as in the interaction of the 
reggeized gluons in the course of the evolution in rapidity. This is rather the initial coupling constant which 
describes the pomeron interaction with the nucleon at zero rapidity. It is unperturbative and not small by 
any estimates. So taking into account higher terms in the expansion of the logarithm in (|19|) is necessary. 
And this is not all. In fact action (I19p and its simplified form (I20p both contain only a single interaction of 
the pomeron with the nucleon. In absence of the pomeron merging and splitting this would give the nucleus- 
nucleus amplitude in the so-called optical approximation, when each nucleon from the projectile interacts 
only with one nucleon from the target and vice versa. Both on physical grounds and inspecting the diagrams 
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without triple pomeron interaction we expect a much richer nucleon-nucleon amplitude, corresponding to 
the full Glauber picture, in which each nucleon from the projectile interacts with every nucleon from the 
target and vice versa. 

So we have to assume that a nucleon from a given nucleus may interact with pomerons many times in 
an eikonal-like fashion. To achieve this the vertex for the pomeron interaction with a single nucleon from 
nucleus A at rapidity y = and transverse point f3, instead of (|15|) . has rather to be 

r^(/3)(e- < ^ t ^ =0 '^ > -l) (21) 

and with with a single nucleon from nucleus B at rapidity Y and transverse point (3 

T B {B - (3) (e-<P<t>iv=Y^)> - i) . (22) 

For the scattering matrix, doing the double sum (I17|) . we obtain, instead of (|18p 

S(Y,B) = (0|( J d 2 /3T A (6)e- < ^ t(s/=0 ' /3)> ) A ( J d 2 /3T B {B - p)e- <p<t>(y=Y ^ > ) B \0), (23) 

which corresponds to a more complicated external action 

A E = AlnJ d 2 [5T A {(i)e- <p ' t> ^ y= ^ )> +B\nJ d 2 /3T B {B - p) e -<P*<3i= Y >P)> . (24) 

If one takes the formal point of view, considers p of to be of the order g 2 and retains only the leading terms 
in g 2 taken to be small, then expansion of the exponentials returns us to the eikonal action (|20p . At small 
rapidities, as desired, the new action (I24p describes the AB amplitude as a full Glauber series and not as 
the optical approximation to it. It is well-known that the optical approximation works not too badly for 
central collisions (B = 0) but quite poorly for more peripheral ones. In QCD the role of peripheral collisions 
is greatly enhanced due to the fact that the non-linear term coming from the pomeron interaction strongly 
damps the contribution of central collisions but does not influence the peripheral ones, which continue to 
grow according to the pure BFKL picture. So the new external action may considerably change the total 
scattering cross-section. As we shall see, performing numerical calculations, the difference introduced by 
transition from (|20p to (|24p is quite large. 

Note that the nontrivial form of the external action (I24p is wholly determined by the nontrivial (3 
dependence of the nuclear distributions in the tranverse plane. For the (unrealistic) case when both Ta and 
T B are constants inside the nuclear volume and B = (central collision) the fields (j) and 4>^ do not depend 
on (3 and one can trivially integrate over (3 in (|24p . Then one returns to the eikonal form of the boundary 
condition ([20]) . 

Also it is worth noting that, for hadron-nucleus scattering, taking into account only fan diagrams 
corresponds to the following procedure. One has to put one of the profile functions, say T B equal to 
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the 5- function, put B = 1, take only the first two terms of the expansion in ft of the logarithm with 
exp(— < pft(y = 0, (3) >) and in (j) and of exp(— < p<p{y = Y, B) >) under the sign of the second logarithm 
and finally take into account only the part proportional to A in obtaining the equation of motion (which 
results in ft = at all values of y). It coresponds to the simple eikonal form (|20l) . with only the term with 
ft remained. 

Differentiation of Ae over the fields leads to the boundary conditions for and ft: 

1 , e -<P* t (s/=0,/3)> 

0(0, k, 0) = —AT A (b)Kp(k) — — . -r-, - -rr (25) 

and 

1 , e -{p<t>(y=Y,0)) 

'' *■ » = 2 BTb(B - mK "^ fW-inrWHM - (26) 



In terms of the gluon densities the external action (|24p acquires the explicit form 



A E = -^Aln U dFhTA^expi — J ^ r(a 2 )ft (0,a 2 M)} \ 

Sin (j d 2 hT B (b - 6l )exp{-^J J ^- T (a 2 )f(Y,a 2 M)}^j , (27) 



iV 2 
2tt 3 



where function r(k 2 ) which characterizes the colour distribution in the nucleon is related to p(k) similarly 
to © : 



r(k 2 ) = ^Vlk'p(k) (28) 



The boundary conditions become 



f(y = 0,fe 2 ,/?) = AT A (P)r(k 2 ) - — — iV - - ^ 2 - ^ — - — (29) 



/ d 2 6' r A (/?') exp{-^ / ^ r(a 2 )/t(0, a 2 . D } 

and 

2tt 3 r da 2 „2\ ffv „2 



e W {-^I^r(a 2 )f(Y,a 2 ,P)} 
J d 2 f3> T B (B - 13') exp{-^ / ^ r(a 2 )f(Y, a 2 , f3<)} 



ft(y = Y,k 2 (3)= BT B (B-P)r(k 2 ) - — - ^ ° ^ q2 ; - — . (30) 



With the fixed Ae, using equations of motion (llip and and a similar equation for one can somewhat 
simplify calculation of the total action, excluding from it, say the free action term, containing derivatives in 
rapidity (see [12]). As a result the total classical action is obtained in the form 

A dass = ~Ai + Ae 

^ [ d *3AT A (8)r(k>) °rt-%SgT(M* 2 M 
N 2 J ' V> f cfp'T A (b')exp{-^f^-r(a 2 )p(0,a 2 ,P')} 
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vr 3 r exp{-^ J d4r(a 2 )f(Y,a 2 ,p)} 

•rW [ d 2 f3 BTb{B — (3) r{k 2 ) , - — ^ M ^ \ >\> - _ . (31) 



N 2 J K; f fPP'TBiB-Wexpl-^fgTWffra?,?)} 

For comparison we note that the eikonal external action ()20|) leads to the boundary conditions which 
are obtained from (|29|) and (|30f) by dropping the exponential and denominator factors. The classical action 
in this case is given just by (l/2)(Ai — Ae)- 



4 Solution 

Calculations of the action and so of the nucleus-nucleus scattering matrix require knowledge of the colour 
dipole distributions p(k) or r{k 2 ) in the nucleon. We take a simple phenomenological ansatz for r 

T ^ = WTW < 32) 

We have reasons to assume that the initial form of the colour distribution of the sources is not so important 
for the behavior of the solution at high rapidities. As it well known from the experience with numerical 
solution of the BK equation, evolution quite fast "forgets" the form of the initial colour distribution. Also, 
as mentioned in the Introduction, in this paper we do not pretend to fit any experimental data. Our aim 
is investigation of general properties of the solution of the nucleus-nucleus scattering amplitude and their 
dependence on the choice of the boundary conditions. Therefore in fact the choice of r(k 2 ) is not so important 
for our calculations. On the contrary, the form of impact parameter distribution of the scattering nuclei is 
very important. Therefore we take the realistic Wood-Saxon parameterization for the profile function: 

Ta{13) = 4^ R\ + a 2 n 2 R A /oo j , ~ 7 -R A +y/W^ s (33) 



where 

R A = 5.7 A 1/3 GeV' 1 , a = 2.725 GeV' 1 . (34) 

In order to numerically solve the equations of motion we used the algorithms described in [12] and in 
|13| . So we shall not repeat technical details of the calculations here, redirecting the interested reader to 
these papers. All calculations, presented below were performed for a s = N c a s /7r = 0.2, with N c = 3 and 
for the case of central collisions (B = 0) of two gold nuclei with A = 197. 

For further comparison we start from a simpler case of 'eikonal' interaction with the external action 
(|20|) . which was previously studied in the papers |12} [T3] but for different choices of initial distributions. 
We recall that in [13] , for scattering of two protons considered as a collection of distributed sources similar 
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to the nucleus, it was found that a symmetry breakdown occurs in the system at sufficiently high rapidity 
Namely at rapidities Y lower than a certain critical rapidity Y c ~ 10 the equations of motion have a unique 
solution which exhibits a natural symmetry between projectile and target: 

f(y,k 2 ) = f\Y-y,k 2 ) (35) 

However at larger rapidities Y > Y c two new solutions are found, which become highly asymmetric as 
rapidity grows: 

fi(y,k 2 )^fl(Y-y,k 2 ), and f 2 (y,k 2 )^ft(Y-y,k 2 ) (36) 
The two solutions are related as 

fi(y,k 2 )=fl(Y-y,k 2 ), (37) 
which guarantees the symmetry between the projectile and target on the whole. 



In the present run of calculations, with the eikonal form of interaction (|20p . nuclei as participants and our 
choice of t(/c 2 ) and Ta(/3), we find the same behavior with practically the same value for the critical rapidity 
Y c ~ 10. At rapidities larger than Y c the minimum of the action is achieved at asymmetric solutions. We 
present one of the asymmetric solution in Fig. for Y=20 and [5 = at maximal rapidities in the evolution 
(y = Y for / and y = for f'). For convenience we plot dimensionless fields h and h' related to / and p 
by 

f(y,k,(3) = ^^k 2 h(y,k,{3), f\y,k,(3) = J^k 2 h\y,k, (3). (38) 

One observes that the larger field h has a form similar to the solution of the BK equation (sum of fan 
diagrams): it has a sharp maximum at large 'saturation' momentum Q s ~ 190 GeV/c. In fact, to a very 
good precision, it is equal to the solution of the BK equation with the same initial condition and so represents 
the sum of only fan diagrams. 

Values of the classical action calculated with the eikonal form of the external action (|20p for both field 
configurations, symmetric and asymmetric, are presented in Fig. 0] as function of rapidity. We see from 
this plot, that at high energy one safely can retain only the non-symmetrical solutions, which provide the 
minimum for the action, and determine the scattering scattering amplitude as 

2 

A el (Y, B = 0) = 1 - ex P ( - Mfi, fl; Y; b}) . (39) 
i=i 

Also one observes that the action rapidly rises up to Y ~ 15 and then saturates at value ~ 1000. Finally we 
note that at high rapidities the values of the action for the asymmetric solutions can be obtained by solving 
first the BK equation for the larger field and then our equation for the smaller field with the found values 
for the larger one. 
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Figure 3: Fields h(y = 20, k, b = 0) and h^(y = 0, k, b = 0) at Y = 20 for the asymmetric solution found 
with eikonal action $2U\) 

Passing to the corrected interaction, with the the external action (|24|) and boundary conditions (|29|) 
and (I30p we first stress that now, from the mathematical point of view, the problem, already non-linear 
and difficult, gets even more complex. The new boundary conditions not only mix values of the fields at 
different points in the participant nuclei, but, which is more important, mix values of the fields at boundary 
rapidities y = and y = Y. So in fact they are not boundary conditions in the proper sense of the word, but 
just complicated non-linear relations between initial and evolved fields at both boundary rapidities. Still 
iterative methods developed in |12| [T3] for the eikonal interaction (j20p proved to be applicable also to this 
new problem. 

The result of our calculations again shows that for rapidities below the critical rapidity Y c the equations 
have a unique solution, which is symmetric under interchange of the projectile and target, that is satisfies 
(|35j) . However the critical rapidity is now much lower: Y c ~ 6. Starting from this value two asymmetric 
solutions appear, which satisfy (|37|) . As compared to the eikonal action, the main difference is that this 
asymmetric pair of solutions is not unique. In fact, as rapidity rises, new asymmetric solutions are found, 
which split from the ones found at lower rapidities. Happily new solutions give greater values for the action. 
The minimal value for it is obtained from the pair of solutions which appears earliest at Y = Y c . The whole 
picture for the ^-dependence of the action obtained from different solutions is shown in Fig. [5j The minimal 
value for the action again first rapidly rises with Y up to Y ~ 6 but then begins to grow very slowly and 
practically saturates starting from Y ~ 10 at A = 270. Comparison with Fig. 0] illustrated in Fig. [6] shows 
that this value is nearly 4 times lower than the one obtained with the eikonal coupling (|20p . The behavior 
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Figure 4: Action obtained with the eikonal external action \2U\) with symmetric (upper curve) and asymmetric 
(lower curve) solutions of the equations of motion. 

of the fields / and f* (or h and h') at high rapidities looks very similar. In Fig. [7] we show the evolved 
fields h(y = Y,k 2 ,f3 = 0) and fv(y = 0,k 2 ,{3 = 0) for one of the asymmetric solutions at Y = 20 (in the 
symmetric case h and h* should be equal) . One observes that field h is very similar to the one obtained from 
the eikonal action. It exhibits a strong maximum at k ~ 270 GeV/c, slightly higher than with the eikonal 
action (c.f Fig. With the growth of Y this maximum shifts towards greater k so that its behavior is again 
similar to the solution of the BK equation, which sums only fan diagrams. As with the eikonal coupling, 
the conjugate field fv is concentrated in the region of small k but it is almost 100 times smaller. Note that 
the contribution to the action comes from the small region of k where the two solutions overlap. Because of 
this the value of the action is very sensitive to the form of the fields. In contrast to the eikonal action case, 
it cannot be reproduced by approximately taking only the sum of purely fan diagrams for one field and the 
expression for the other resulting from its evolution (which would return us to the eikonal action case). 



13 



Action 

1600 - 



400 



Au-Au,B=0 



1200 



I 000 



800 



600 



400 



200 




Figure 5: Action obtained with new boundary conditions and \29\) . The uppermost curve shows the 
action for the symmetrical solution. Lower curves correspond to the family of asymmetrical solutions. The 
upper curve from this family, which steadily grows with Y coincides with the one obtained from the old 
eikonal action jl20\) (c.f Fig.\^. 

5 Discussion 

In this paper we applied the effective field theory of interacting BFKL pomerons to nucleus-nucleus scat- 
tering. Our calculations were performed with two different boundary conditions. The difference in them is 
that the new, corrected boundary conditions take into account additional diagrams, physically relevant but 
missed with the simpler 'eikonal' variant of boundary conditions. The new initial conditions are nonlocal on 
the nuclear scale, non-linear in fields and mix the initial and fully evolved fields. They allow for the eikonal 
type interaction of each nucleon from the nuclei with pomerons. 

Numerical calculations with both types of boundary conditions revealed a certain similarity but also 
definite difference. The similarity is in the general structure of the obtained solutions of equations of 
motion. In both cases there are symmetric and asymmetric solutions, and at high energies only asymmetric 
solutions need to be taken into account. The difference is that in the case of the new boundary conditions 
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Figure 6: Action obtained for the leading asymmetrical solutions of the equations of motion for the old (upper 
curve) and new (lower curve) initial conditions. 

not one, but a whole family of asymmetric solutions is found. Fortunately, of all these solutions only one 
is leading and provides the minimum of the action. All the rest give much greater values of the action and 
therefore can be neglected. Comparing this minimal value with the one obtained with the eikonal initial 
conditions, one observes that additional diagrams like shown in Fig. [8^,, which are included by the corrected 
boundary conditions, provide a much smaller minimum. They make the nuclei substantially greyer in the 
high-energy limit. Note that in the description of DIS by the BK equation, the need for initial conditions 
which include multiple interaction of each nucleon was stressed already in |10| . Applied to hadron- nucleus 
scattering, initial conditions of [10j generate diagrams crudely described by Fig. [8b (in fact the lowest 
lines should be not reggeized but simple gluons). In our approach hadron- nucleus scattering also involves 
diagrams like Fig. [5^, (for DIS one can prove that such diagrams are in fact included into Fig. [5b). 

We have found that both for the old and new boundary conditions one of the pomeron fields behaves 
very similar to the solution of the BK equation while the other is very small. This of course raises a question 
whether in our nucleus-nucleus case it is possible to retain only fan diagrams and thus radically simplify 
the problem reducing it to the BK equation for only one field. Unfortunately this is not possible, at least 
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Figure 7: Fields h(y = 20, k, b = 0) and h){y = 0, k, b = 0) at Y = 20 /or i/ie solution which minimizes the 
action. 

for the nucleus-nucleus amplitude, since, as mentioned, the contribution in fact comes from a small region 
where the two fields overlap. However this structure of the solution may simplify the study of the double 
inclusive cross-section, known to be given by a host of complicated diagrams (in contrast to the single 
inclusive cross-section). The fact that the smaller field is concentrated in the region of very small momenta 
may imply that the bulk of the contribution at physical momenta comes from the greater field and thus 
from pure fans, tractable via the BK equation. We leave this problem for later studies. 

Lately there have been some progress in introducing the running coupling into the BFKL scheme |20ll21l 
[2"2"| I23). In [21] (and earlier in [55], where the running coupling was introduced phenomenologically) it was 
found that the running of the coupling significantly lowers the rate of growth of the saturation momentum 
with rapidity, thus bringing it into agreement with experimental findings. In [23] the form of the triple 
pomeron vertex with the running coupling has been found. So, in principle, the door is open to generalization 
of the effective pomeron theory and its application to nucleus-nucleus scattering with the running coupling. 
However the structure of the found vertex is most complicated and the equations which are expected for 
the pomeron fields do not seem easily tractable. Still we do not expect much change in the overall behavior 
of the solutions, except that both the action and saturation momentum for the larger field are expected to 
become lower. 

Finally a couple of words about the loops. In the nucleus-nucleus case they are suppressed by factors 
A 1 ! 3 and B 1 / 3 coming from the interaction with nucleons. They are not suppressed at all in purely hadronic 
interactions and there the problem is acute. A crude estimate of a single pomeronic loop contribution was 
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Figure 8: The example of diagrams missed with the eikonal form of action and included in calculations with 
more sophisticated boundary conditions. 

made in [26J, with the conclusion that it was negligible at reasonable rapidities. More instructive results 
have been achieved only for the zero-dimension variant of the local pomeron RGFT [14} [T5| [To] \13\ \T7\ [28] 
Surprisingly, the obtained full quantum solution for the amplitude of two interacting one dimensional "pro- 
tons" turned out to be close to the amplitude for the asymmetrical solution of the equations of motion. 
We do not know, if such a property will be preserved in the theory of interacting QCD pomerons. Up 
to now there are no full QCD loops calculations similar to the ones performed for zero dimensional mod- 
els. In principle, there exists a complete and consistent theoretical framework for loop calculations in the 
nucleus- nucleus case ( |18j). However its practical realization seems to be only possible by numerical lattice 
calculations, which lie beyond our present interests. 
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